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ABSTRACT 

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin 
in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a 
spin change can be characterized by a stochastic random walk of angular momentum. Also, spin 
distributions generated by random walks are fairly consistent with those directly obtained from N- 
body simulations. We derived a stochastic differential equation from a widely used spin definition and 
measured the probability distributions of the derived angular momentum change from a massive set of 
halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving 
spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are 
found to influence the angular momentum change. The spin distributions generated in the mean-field 
or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, 
which seems to be caused by the correlated random walks. We verified the assumption of randomness 
in the angular momentum change observed in the N-body simulation and detected several degrees of 
correlation between walks, which may provide a clue for the discrepancies between the simulated and 
generated spin distributions in the voids. However, the generated spin distributions in the group and 
cluster regions successfully match the simulated spin distribution. We also demonstrated that the 
log-normality of the spin distribution is a natural consequence of the stochastic differential equation 
of the halo spin, which is well described by the Geometric Brownian Motion model. 

Subject headings: cosmology: theory — halo spin — large-scale structure of universe — methods: 
numerical — methods: statistical — random walk 


1. INTRODUCTION 

Primordial angular momentum is believed to de¬ 
cay with the expansion of the universe and the ob¬ 
served rotations of galaxies and clusters are believed 
to originate from tidal interactions with local environ¬ 
ments (the tidal torque theory) or the acquisition of 
orbital angular momentum from infalling matte r (the 
merging theory). In the tidal torque theory (iHovlel 

I94ft[^eblesll96^ll^e fc Penl200il : lProciani et alJl2002t 

Lee fc Erdogdd 120071 1. a collapsing cloud at the turn¬ 

around epoch may obtain a tidal torque due to the 
misalignment between the local tidal shear and struc¬ 
tural inertia tensor. Consequently, the spin directions 
of halos may provide information on the nearby mass 
distri b ution (I Ara gon-Cal vo fc Yan^ I20l4 iZhang et al.l 
I2009I iTrowland et al.l 201,111. although some studies 
l|Lee fc Erdogdull2007t iVarela et al.ll2012ll have reported 
inconclusive or opposing results from o bservations. On 
the other hand, the merging theory ([Gardner] 120011 : 
IVitvitska et al.l l2002t iMaller et al.l 1200211 focuses on the 
change in angular momentum in hierarchical clusterings 
through the off-center infall of nearby o bjects/matter 
to dr ive the sudde n change in halo spiii (iBett &: FrenkI 
1201211 . However, iD’Onghia fc Navarrol (|2007ll argued 
that simulated halos do not show any obvious correla¬ 
tion between spin and merger history. 

Each theory views the same object from different an¬ 
gles. A dark matter halo is not static but dynamically 
evolves through the sometimes violent, sometimes quiet 
ingestion o f nearby matter. According to hiera rchical 
clustering (jCole fc lacevi Il996t I Jiang et^ 1201411 , large 


halos grow in mass by assimilating nearby smaller halos. 
However, a small satellite halo tends to have a spiral 
infall motion with a finite orbital angular momentum 
around the central halo, which, after merging, causes 
a net change in rotational angular momentum of the 
merging remnants. This non-radial infall motion of satel¬ 
lites may be derived from the tidal torque of anisotropic 
matter distributions around the central halo. However, 
the acquisition of satellite orbital angular momentum de¬ 
scribed in the merging theory is analogous to that of a 
mass shell in the tidal torque theory because they are 
based on the same ground: the tidal torque generated by 
local environments. The only difference is whether they 
consider the infalling source to be a satellite or a mass 
shell. By extending the original tidal torque theory of a 
spherical collapse to the multiple-shell collapses at dif¬ 
ferent epochs, the tension between the two theories may 
be relaxed. 

The distributions of halo spin have been reported 
to be well fitted by the log-normal function (iGardnerl 
20011: iBailin fc Steimnetzl 120051 iBullock et al. I 120011 : 
Muhoz-Giiartas et al.l I20il1 l. with small dev iations de¬ 
pend i ng on the definition of the halo samples (IBett et al. 
2007 : iKnebe fc Power! 120081 lAntonuccio-Delogu et al. 
201011 . Researchers have found that massive cluster halos 
are not yet relaxed from recent mergers and are responsi¬ 
ble for the over-population o f the high-A tail compared to 
the lo g-normal distribution (jBett et al.ll2007t iHahn et al.l 
[ 21 ^. _ 

lAntomiccio-Delogu et al.l 1)2010( 1 explored the origin 
of the log-normality of spin distributions based on the 
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stochastic variables of halo angular momentum (J), total 
energy (E), and mass (M). Theirs was the first attempt 
of its kind to address this issue. They introduced stochas- 
ticity between the two terms | and E and argued 

that deviations from the log-normal functional form are 
mainly due to correlations between the two terms. How¬ 
ever, they did not provide a physical or probable reason 
for the stochasticity between the two terms. 

In the present study, we analyze the origin of the 
halo spin distribution. Using simulated mass histories 
of individual halos and the stochastic differential equa¬ 
tion of the angular momentum change, we model the 
spin evolution using the stochastic models distinguish¬ 
ing between the roles of major merging and accretion 
i n shaping the spin di s tribut i ons in various situations 
(iHetznecker &: BurkertI 120061: lAntonuccio-Delogu et al.l 

I2nn2t inaTdneHl2nniU T, emson fc KauffmannI 1 19991 1. 

This paper is organized as follows. In section [H we 
describe the stochastic for spin evolution. Then, in sec¬ 
tion [21 we describe the simulation data, merging trees, 
and local environments of a halo. From the simulated 
merging history, we measure distributions of angular mo¬ 
mentum change with respect to halo mass, spin, local en¬ 
vironments, infall mass ratio, and redshift in section 0) 
In section [SI using the simulated mass histories of in¬ 
dividual halos, we randomly generate spin changes and 
compare the spin distributions of the stochastic model 
and 7V-body simulations. In section |6l we verify how a 
random-walk process produces the log-normal distribu¬ 
tion of the halo spin. We summarize the results of our 
study in section [7] 

All of the logarithms used in the equations and figures 
in this study are base 10. 


may show a signihcant discrepancy between A and Aq be¬ 
cause most are not fully virialzed and their radial profiles 
do not follow the isothermal model. 

The gravitational potential of a halo should be mea¬ 
sured using the same finite softening length (e) as 
adopted in iV-body simulations because the chosen po¬ 
tential estimation strongly depends on the potential 
model: the Newtonian (1/r) or Plummer (l/Ve^ -I- r^) 
model. The d ifference increases at higher redshift 
(|Ahn et al.l l2014f) . To comply with the simulated po¬ 
tential, we measure the halo potential using 


<I>(r) 


GmiTTij 


( 4 ) 


where e = 36 /i“^kpc. Therefore, great care must to be 
taken when interpreting the halo spin value. 

In this study, we adopt a new spin parameter A' 
(lAhn et al.ll20I4il . 


A' = /WA' 

/ 0.26 


( 5 ) 

( 6 ) 


whi ch is introduced t o correct for all possible systematics 
(see lAhn et ahl 120141 . We call f{z) a correcting factor. 
Hereafter, whenever we refer to a halo spin, we mean the 
corrected spin and use A to denote it. Of course, one may 
return to Ag simply by dividing the spin value by f{z). 


2.2. Stochastic Model: Random Walk of Halo Spin 

From equation (l2|), we derive a stochastic differential 
equation of spin change as 


2. ANALYTIC DESCRIPTIONS 
2 . 1 . Definition of Spins 


iPeeblesI (|I969ll first proposed a dimensionless param¬ 
eter to quantify the rotation of a cosmic object or halo 
as 


GM5/2 ’ 


( 1 ) 


where A is the spin, E is the total energy, M is the virial 
mass, and J is the rotational angular momentum of a 
halo. However, measuring halo potential is practically 
bothersome, and the potential of an open system is some¬ 
times ill- defined. Therefore, an other definition was pro¬ 
posed bv iBullock et al. I (1200111 . which reads as 


A'o = -7^ -, (2) 

V2MRV 

where R is the virial radius, and V is the virial circu¬ 
lar velocity. The virial radius can be derived from the 
definition of virial mass as 

47r o 

M = —R^pc{z)A^ifiz), (3) 

where pc is the critical density at redshift z, Avu — 
(IStt^ + 82a: — 39x^1, an d x = fl{z) — 1 for a flat uni¬ 
verse (jBrvan fc Norman! 1199811 . These widely used spin 
parameters are equivalent to A in the case of an idealized 
halo; halos are fully virialized and their density profiles 
follow the isothermal model. However, simulated halos 


c? logic A ^ cl logic \J\ , / 1 5 

dlogiQM dlogi^M ^ 3’ 

where 

^ 1 cHogio(pe(^:)Avir) dlogic/(^^) 

6 dlog^gM dlog^gM 


( 7 ) 

( 8 ) 


Here, a{z) is a cosmological effect depending on the red¬ 
shift and model parameters. However, its contribution 
is relatively trivial compared with the constant factor of 
—5/3 in equation d?]) unless Idlog^g M\ < 10“^. This im¬ 
plies that the background history of the universe does not 
much affect the evolution of the spin parameter (Lemson 
& Kauffmann 1999; Cervantes-Sodi et al. 2008). 

This differential equation was derived to model the spin 
change as a function of the mass event (mass-merging 
or mass-loss event) accompanying a change in angular 
momentum amplitude | J|. Here, we introduce a Markov 
chain process of spin evolution as 


A logio A 

A logio Ad 


D{M, A, 2 , etc.) — Dc, 


( 9 ) 


where Dc = 5/3 — a{z) and the angular momentum 
change rate is defined as 


D{M, A, z, etc.) 


A logio 1^1 

A logio Ad ■ 


( 10 ) 


Since (x{z) is relatively small, we will not explicitly men¬ 
tion it, even though we include this contribution in the 





































3 


TABLE 1 TABLE 2 

Simulation characteristics Basic information on the merging data 


Np 

-^box 

A^step 


fib 

Da 

h'° bs^ 

2048^ 

737.24 

120 3000 

0.26 

0.044 

0.74 

0.72 1.26 


® Starting redshift of the simulation 

^ Hubble expansion parameter divided by 100 kms“^Mpc“^ 
Bias factor at i? = 8 /i“^Mpc 


subsequent calculations. We stress the importance of 
a{z) later in Appendix lAl 

Throughout this stochastic approach, our primary as¬ 
sumption is that a spin change, Alogj^pA, does not de¬ 
pend on the past history of a halo but can be completely 
described by the current physical status of the halo; the 
spin change is determined by the stochastic probability 
of angular momentum change, P{D), which will later be 
shown to be a function of M, A, Alog^gAf, z, and lo¬ 
cal environments. Later, we will prove that if P{D) is 
a Gaussian distribution, then the resulting spin distribu¬ 
tion follows the log-normal shape. In the stochastic equa¬ 
tion, there is one independent variable, A logj^g M, which 
should be determined outside of the model. The halo 
mass history could be provided by A-body simulations 
or by an extended Press & Schechter (EPS) formalism. 
However, in this analysis, we limit our study to the simu¬ 
lated mass histories. Hereafter, we define a major merger 
event and an accretion (or minor-merger) event as mass 
changes of Alog^oAf > 0.1 and 0 < Alog^gM < 0.1, 
respectively. 

3. DATA 
3.1. Simulations 

We run a series of cosmological A-body simula- 
tions with 2048^ pa rticles using the GOTPM code 
(|Dubinski et al.l [20041 1. Four of the simulations are ba¬ 
sically identical but were run with different sets of ini¬ 
tial conditions. The simulations have a cubic simulation 
volume with a side length of Tbox = 737.28 /i“^Mpc. 
The cosmological model adopted in the simulations is a 
WMAP 5-year cosmology. The fractions of matter, bary- 
onic matter, and dark energy are 0.26, 0.044, and 0.74, 
respectively. The initial power spectrum is obtained from 
the GAME Source packageQ- Table [T] summarizes several 
characteristics of the simulations used in this paper. 

We applied the Friend-of-Friend (FoF) method to iden¬ 
tify groups of particles with a percolation linking length, 
^FoF, which is usually (also in this paper) set to be 0.2 
times the mean particle separation. This ensures that 
the identified halos to have a mean over-density of 180, 
satisfying the cosmic virialization conditions in the spher¬ 
ical top-hat collapse model. The minimum mass of a halo 
identihed with 30 member particles is about 10^^ 

Data on the member particles of each halo are stored at 
44 redshifts during the simulation run. Table [2] contains 
some basic information on the merging data extracted 
from these A-body simulations, including the simulation 
step number {first column), redshift {second), time be¬ 
tween merging steps {third), number of halos {fourth), 
number of seed halos {fifth), and mean halo number den¬ 
sity {last). Here, seed halos are introduced to define a 

^ http://camb.info/sources 


step^ 

z 

AA 

-^samp^ 

A^bkg'^ 

n** 

250 

10 

- 

174 

174 

4.3 X 10“’^ 

278 

9 

125.7895 

7,317 

713 

4.6 X 10"® 

353 

7 

336.9363 

264,329 

47,823 

0.00016 

407 

6 

242.5939 

1,135,317 

3,735 

0.00071 

440 

5.5 

148.2521 

2,142,976 

11,608 

0.0013 

479 

5 

175.2068 

3,841,831 

34,147 

0.0024 

501 

4.8 

98.8345 

5,025,833 

56,128 

0.0031 

525 

4.5 

107.8197 

6,471,073 

89,473 

0.004 

551 

4.3 

116.8045 

8,190,543 

139,552 

0.0051 

580 

4 

130.2820 

10,239,195 

213,106 

0.0064 

605 

3.8 

112.3119 

12,072,304 

292,010 

0.0075 

633 

3.6 

125.7898 

14,159,663 

397,941 

0.0088 

662 

3.4 

130.2818 

16,313,976 

526,000 

0.01 

695 

3.2 

148.2521 

18,714,095 

692,528 

0.012 

731 

3 

161.7294 

21,221,990 

893,728 

0.013 

771 

2.8 

179.6993 

23,841,886 

1,137,628 

0.015 

815 

2.6 

197.6691 

26,477,254 

1,419,370 

0.017 

865 

2.4 

224.6239 

29,135,986 

1,746,844 

0.018 

920 

2.2 

247.0869 

31,672,320 

2,105,555 

0.02 

983 

2 

283.0264 

34,107,018 

2,500,701 

0.021 

1018 

1.9 

157.2365 

35,262,139 

2,708,351 

0.022 

1055 

1.8 

166.2222 

36,346,533 

2,918,928 

0.023 

1095 

1.7 

179.6993 

37,374,958 

3,132,572 

0.023 

1138 

1.6 

193.1770 

38,330,932 

3,347,949 

0.024 

1185 

1.5 

211.1466 

39,224,698 

3,565,667 

0.024 

1235 

1.4 

224.6238 

40,019,836 

3,777,393 

0.025 

1290 

1.3 

247.0869 

40,747,200 

3,985,683 

0.025 

1350 

1.2 

269.5488 

41,384,587 

4,188,379 

0.026 

1415 

1.1 

292.0116 

41,921,471 

4,379,074 

0.026 

1487 

1 

323.4587 

42,370,158 

4,561,970 

0.026 

1567 

0.9 

359.3986 

42,728,496 

4,731,685 

0.027 

1656 

0.8 

399.8306 

42,996,082 

4,887,605 

0.027 

1754 

0.7 

440.2634 

43,161,903 

5,027,849 

0.027 

1866 

0.6 

503.1586 

43,252,426 

5,155,748 

0.027 

1992 

0.5 

566.0531 

43,252,029 

5,264,860 

0.027 

2136 

0.4 

646.9167 

43,180,668 

5,362,503 

0.027 

2302 

0.3 

745.7516 

43,035,938 

5,443,524 

0.027 

2496 

0.2 

871.5429 

42,835,401 

5,513,029 

0.027 

2725 

0.1 

1028.778 

42,586,126 

5,570,734 

0.027 

2776 

0.08 

229.1170 

42,533,514 

5,581,988 

0.027 

2856 

0.05 

359.3983 

42,446,174 

5,597,921 

0.026 

2941 

0.02 

381.8615 

42,359,628 

5,611,664 

0.026 

3000 

0 

268.061 

42,298,156 

5,621,397 

0.026 


Simulation time step in units of Myr 
^ Time interval between steps 

Total number of halos with members equal to or larger than 30 
simulation particles 

Number of background halos needed to trace density field 
® Mean number density of halos in units of /i®Mpc“® 

local density which will be described in detail in sec¬ 
tion [331 

3.2. Main Merging Tree 

To build halo merging trees, we used information for 
member particles of each FoF halo. If a “majority” of 
member particles of a halo are also members of a halo at 
the previous time step, then we link them with a mainr- 
merging tree line. We call the former halo an “ancestor” 
and the latter a “major descendent.” If a halo has no ma¬ 
jor descendent, then we terminate the line. The chains of 
these merging lines consist of a merging tree with which 
we are able to trace the histories of halo angular momen¬ 
tum, mass, and other physical quantities. One should 
note that the main merger tree is a trimmed version of 
the merging tree widely quoted in other studies. 

Note that, on average, about 25% of halos experience a 
mass loss in a time step of this study, which means that 
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Fig. 1.— Several examples of spin evolution of simulated halos. 
Each color represents the spin trajectory of a single main-merging 
tree. For the x axis, we use halo mass rather than time or redshift. 
Halos at 2 = 0 are chosen with mass 10^® < M < 2 x 

its major descendent (or successor) becomes less massive. 
This happens when a halo undergoes a violent merging 
process and is broken into irrultiple less-massive remnants 
or when a halo passes another and is stripped of some 
member particles. 

3.3. Spin Evolution of Simulated Halos 

Now, let us look at the spin evolution of simulated 
halos and their main-merging trees. In Figure [U we 
show several typical examples selected from the simu¬ 
lations. There are several significant spin decreases to 
less than A = 0.01 or increases to greater than 0.07, but 
halo spins immediately bounce back to settle over the 
range of 0.01 ^ A ;< 0.07. Here, we suppose that spin 
change is random, but there is a hidden mechanism that 
regulates the spin value to lie in a certain range. 

A sudden massive merger event may cause a halo to ro¬ 
tate quickly, resulting in a substantial loss of angular mo¬ 
mentum due to the ejection of high-angular momentum 
material through violent relaxation. A smaller halo may 
pass a nearby larger halo (being attached when they en¬ 
counter but separated afterward), after which point the 
larger halo experiences a sudden increase and decrease in 
spin value. Of course, during the encounter, there could 
be an exchange of finite angular momentum and the host 
may not fully recover its original angular momentum. As 
stated before, several events of mass loss events due to 
the separation of satellites can be seen in Figure[T] (a line 
segment moving to the left). 

3.4. Local Environment 

We define a local density measured from 10 nearest 
neighbor halos as a halo environment. The halos tracing 
the local density, called “background halos”, are all halos 
with mass greater than Mg = 10 ^^ (equal to 

the sum of 300 particle masses) at z < 7. However, at 
higher redshifts (z > 7), the number of massive halos is 
much lower and, therefore, we reduce the mass limit to 
Mg = 2 X lO^^/i“^M 0 at 7 < z < 9 and even further to 
Mg = 10^^/i“^Mq at z > 9. 


The variable-width spline kernel is adopted to adap¬ 
tively resolve high-density regions as 

Ar„=io 

E (11) 

i=i 

where is the distance from the ith halo to the jth 
background halo, and E is the distance to the tenth 
nearest neighbor. The density measure was first used 
with 20 neighbors in Park et al. (2007). The kernel size 
varying with environments helps dense regions not to be 
oversmoothed. 

The spline kernel has the form 

f 1 — 6x'^ + if 0 < a: < 

W{x) = — < 2(1 — x)^ if ^ < X < 1, and 

^ [ 0 if X > 1, 

where x = r/h. The median radius of the kernel is h = 

7.5 h-lMpc. 

The local environment is parameterized by 

Apio(r-) = (12) 

P 

where p is the mean density of the background ha¬ 
los, which is measured to be 3.54 x 10“^(h“^Mpc)“^ = 
(6.56h“^Mpc)“^. If a region has a density of Apio < 0.7, 
then we refer to it as an under-dense or void region. Re¬ 
gions 0.7 < Apio < 2, 2 < Apio < 10, 10 < Apio < 100, 
and Apio > 100 are referred to as mean-field, group, 
cluster, and highly clustered regions, respectively. 


4. PROBABILITY DISTRIBUTION OF ANGULAR 
MOMENTUM CHANGE, P(D) 

4.1. General Descriptions of P{D) 

We measured the distribution of angular momentum 
changes, P{D), as a function of M, A, Alog^gM, z, 
and Apig. First, we divided all of the merging events 
into five-dimensional subsamples stretching over the en¬ 
tire parameter space, and measured the distribution of 
D in each subsample. The D was derived by estimat¬ 
ing the changes in mass and the angular momentum 
(Eq. (To]). Between merging steps (ti and ^ 2 ), the log 
mass change would be log^g M 2 — log^g Mi and the an¬ 
gular momentum change in the log scale is obtained as 
logio 1 ^ 2 ] — logio where subscripts 1 and 2 indicate 
the time step. Then, we fit the distribution with a bi- 
modal Gaussian function as 


P{D-,p,i,ai,p2:Cr2) 


h 

\/27rcrf 


exp 


1-/1 


: exp 


1 {D-P2f 

2 


2 a\ 


(13) 


where p is the Gaussian mean, and <j is the standard de¬ 
viation. The overall mean {po) and standard deviation 
{(Jd) are obtained as 

pD = fipi + — fl)p2 (14) 

crD=fi {erf+ dl) + {1-fi) (aj + df) , (15) 

where di = pi — po and d2 = P 2 — pD- 
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TABLE 3 

Parameter, ranges'* applied to the measurement of 
P{D) 


id 

halo mass (M 12 ) A 

Alog^g M 

Apio 

Z 

1 

0.1 

0.004 

-0.3 

0.7 

5 

2 

0.3 

0.006 

-0.1 

2 

4 

3 

0.5 

0.0085 

-0.07 

10 

3 

4 

0.8 

0.01 

-0.05 

100 

2 

5 

1 

0.015 

-0.03 

00 

1.5 

6 

2 

0.02 

-0.02 

- 

1 

7 

4 

0.025 

-0.01 

- 

0.6 

8 

5 

0.03 

-0.0065 

- 

0.4 

9 

7 

0.035 

-0.004 

- 

0.2 

10 

10 

0.038 

0 

- 

0.0 

11 

20 

0.04 

0.004 

- 

- 

12 

50 

0.043 

0.0065 

- 

- 

13 

500 

0.045 

0.008 

- 

- 

14 

5000 

0.047 

0.01 

- 

- 

15 

00 

0.05 

0.015 

- 

- 

16 

- 

0.053 

0.02 

- 

- 

17 

- 

0.056 

0.03 

- 

- 

18 

- 

0.06 

0.05 

- 

- 

19 

- 

0.063 

0.08 

- 

- 

20 

- 

0.067 

0.1 

- 

- 

21 

- 

0.07 

0.15 

- 

- 

22 

- 

0.072 

0.2 

- 

- 

23 

- 

0.075 

0.3 

- 

- 

24 

- 

0.078 

0.7 

- 

- 

25 

- 

0.08 

00 

- 

- 

26 

- 

0.09 

- 

- 

- 

27 

- 

0.1 

- 

- 

- 

28 

- 

0.11 

- 

- 

- 

29 

- 

0.12 

- 

- 

- 

30 

- 

0.13 

- 

- 

- 

31 

- 

0.15 

- 

- 

- 

32 

- 

0.17 

- 

- 

- 

33 

- 

0.2 

- 

- 

- 

34 

- 

0.25 

- 

- 

- 

35 

- 

0.3 

- 

- 

- 

36 

- 

0.4 

- 

- 

- 

37 

- 

0.5 

- 

- 

- 

38 

- 

0.7 

- 

- 

- 

39 

- 

1.2 

- 

- 

- 

40 

- 

00 

- 

- 

- 


^ The table lists the upper bound of the parameter range 
except the redshift where we instead write the lower bound. 
Then, a parameter range is pi—\ < P < Pi-, where P is a 
parameter, and pi is a value written in the ith row. For po? 
the lowest possible value is assumed. 

Figure [2] shows an example of the probability distri¬ 
bution of angular momentum change in a void sample 
(Apio < 0.7) with parameter ranges 0.03 < A < 0.035, 
0.01 < AlogipM < 0.015, and 1 < M 12 < 2, where 
Mi 2 = M/(10^^ a fitting to the bimodal Gaus¬ 

sian is shown by a solid line. The star symbol with a 
solid error bar in the middle of the plot represents the 
mean with l-cr standard deviation, while the lower two 
error bars (dotted and dashed) mark those of each Gaus¬ 
sian component. Also, we overplot the single Gaussian 
fit (long dashed curve), which shows a poorer fit to the 
fat tails on both sides (more easily seen in the upper 
panel). Spin distributions generated based on the sin¬ 
gle Gaussian fitting function show large deviations from 
simulated distributions, while the bimodal Gaussian has 
a better result in describing simulated spin distribution 
at most redshifts. 

There are several points to note in this plot. First, 
the distribution around a peak is well described by a 
single major Gaussian, while the fat tails on both sides 
require an additional minor broad Gaussian. Second, the 



0.035, 0.01 < Alog^^pM < 0.015, and Apio < 0.7 in the redshift 
interval of 0 < z < 0.2. A total of 2745 mass events are recorded in 
these parameter spaces. Bottom: The solid line with symbol and 
error bar shows the least fit, mean, and standard deviation of 
the distribution. Also, we plot the mean and standard deviation 
of each Gaussian component using dotted (major component) and 
dashed (minor component) lines. Their corresponding Gaussian 
components are shown using the same line styles. The long dashed 
line is a single Gaussian fit to the symbols. Top: the deviations 
of P{D) with respect to the bimodal Gaussian fit are shown with 
circles. 

peak positions of the Gaussian components are not well 
separated from each other, implying that the combined 
distribution is somewhat symmetric. Most subsamples 
show this kind of distribution, although several show 
significant deviations from a symmetrical shape. Such 
asymmetric distributions of D may play a crucial role 
in distorted cluster spin distributions from a log-normal 
form, which will be discussed in subsequent sections. 

4.2. Comparison of P{D) among Halo Subsamples 

In Figure [3l we show the redshift dependence of P{D) 
in halo subsamples for a small mass infall ratio (0.01 < 
Alog^o Af < 0.015) and an average spin value (0.038 < 
A < 0.04) in the mean field (0.7 < Apio < 2). As the 
redshift approaches zero, the standard deviation of the 
distribution becomes narrower and the mean value be¬ 
comes larger. From this behavior, we understand that 
the infall direction is more likely to be random at higher 
redshift because /iu ~ 0 (see the distribution color-coded 
in blue at redshift between z = 3 and 4). At a lower red¬ 
shift, angular momentum changes tend to be larger, sug¬ 
gesting that the orbital angular momentum of infalling 
matter tends to be a bit more aligned with (or positively 
biased to) the halo rotational axis (also see Appendix 
[0. If 11 ~ Dc {vertical dotted line), then we expect that 
there would be no substantial change in the spin value. 

Figure m shows a similar plot, but this time with slowly 
rotating halo samples in the accretion mode. All dis¬ 
tributions have pd > Dc indicating that low-spin (or 
slow-rotating) halos are highly likely to increase in spin 
value. The distribution spread is larger than in the mild- 
rotating samples (see Fig. [3]), implying that the spin 
change is larger for slow-rotating halos. This is because 












6 



-25 -20 -15 -10 -5 0 5 10 15 20 25 

DsdlnlJj/dlnM 


Fig. 3.— Redshift evolution of P{D) for a spin range of 0.038 < 
A < 0.04 (in a moderately rotating case). Other parameter ranges 
of the sample are shown in the legend. The vertical dotted line 
marks D = Dc as a reference. 



Fig. 4.— Same as Fig. [3] but for low-spin sample of 0.015 < 
A < 0.02, which is used to isolate the spin effect on the angular 
momentum change. 

incoming orbital angular momentum may overwhelm the 
previous small rotational angular momentum. 

Next, we investigate the effect of environment on D. 
Figure [5] shows that halos in cluster environments experi¬ 
ence stronger angular momentum changes inferred from 
larger tails (or higher probability of large-Z? changes). 
This means that the angular momentum changes are 
statistically greater than those in lower-density regions, 
or that biased (or systematic) accretion or high-velocity 
mass infall is more frequent in cluster regions. Also, one 
should note that each Gaussian component has a similar 
peak position in most of the fittings. However, when we 
select a sample of higher local density, the centers of two 
Gaussian components begin to separate from each other. 
This effect on spin distribution will be described in detail 
in section and 15.31 

Figure [S] provides an overall view of the effect of mass 
growth (Alogj^gM) on angular momentum change in 
samples of average spin values (0.035 < A < 0.038). The 
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Fig. 5.— Environmental dependence of P(D) for the moderately 
rotating halo samples used in Fig. [3] 
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Fig. 6. — Dependence of P(D) on the amount of mass infall ratio 
(Alogj^Q M) using samples of 1 < M 12 < 2, 0.7 < Apio < 2, and 
0.035 < A < 0.038. 

average value of D (or ^d) does not deviate much from 
Dc (vertical dotted line) except in major-merging cases 
(0.15 < Alogj^gM < 0.2). The major merger signif¬ 
icantly increases halo spin, while accretion (or minor- 
merger) is unlikely to produce a significant systematic 
change in spin for a halo. The standard deviation in 
the accretion mode seems to be anti-correlated with 
A logj^g M, but the significance of the anticorrelation is 
a bit relaxed when weighted with the mass change (like 
D X Alog^o AZ) to update the halo spin (see Eq. [9]). 

4.3. Asymptotic Walk of Halo Spin 

More interesting features can be obtained if one inves¬ 
tigates the dependence oi piu on sample spin. Figure [7] is 
a typical example showing how pLjo changes with sample 
mass and spin. By changing the sample spin while fixing 
other parameters, we determine the asymptotic state of 
the spin distribution, li fijj > Dc, then the halo spin will 
increase; however, if fijo < Dc, then the spin tends to de¬ 
crease. If /iD has a negative slope and crosses Dc, then 
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A 

Fig. 7. — Changes oi as a function of sample mass in case 
of accretion (0.01 < Alogj^gM < 0.015). The horizontal line is 

fin = Dc- 

the spill random walk (Eq. [5|) may converge to a charac¬ 
teristic spin value where crosses D^- If we define Ac as 
a characteristic spin as /i_d(Ac) = Dc, an the asymptotic 
state of the spin distribution is produced with a width 
depending on the shape of P{D) at A = Ac. Virtually 
every mass sample shows a negative slope, and the range 
of the characteristic spin is 0.01 ^ Ac < 0.05. A more 
massive sample tends to have a lower value of Ac- 

If /Id <0, then the rotational axis (•/) of the halo and 
the direction of the orbital angular momentum (L) of 
infalling matter tend to be anti-correlated. However, if 
0 < < Dc, then those two vectors show a positive cor¬ 

relation, while the mass increase is more dominant than 
the angular momentum increase; as a result, the halo 
spin decreases (see Eq. HD. li fiD > Dc, then these two 
vectors show a strong correlation, and halo spin tends 
to increase accordingly. Therefore, a positive correlation 
between two angular momentums does not always in¬ 
crease the halo spin depending on the amount of angular 
momentum change. 

The behavior of hdW may shape a so-called asymp¬ 
totic state of the spin distribution. If a halo spin is sig- 
nihcantly larger or smaller than Ac, then the halo may 
experience a kind of pressure toward Ac; as a result, the 
halo spin tends to move toward Ac. Therefore, the value 
of Ac plays an important role in guiding the spin evolu¬ 
tion. Of course, all halos do not immediately jump to 
Ac but show a spread over a certain range of A around 
Ac, which is regulated by the shape of P{D\X = Ac). 
More massive halos tend to have lower Ac, which means 
that lower spin values are preferred in more massive halo 
samples. This is consistent wit h recent findings for mass- 
dependent spin distributions (lAntonuccio-Delogu et al.l 
[2nint iKnebe fc Poweill200a iBett et al.ll2007l) . 

EigurelH demonstrates how the characteristic spin (Ac) 
has varied with the infall mass ratio (A logj^Q M) in recent 
epochs (0 < 2 < 0.2). Broadly speaking, the halo spin 
tends to increase as the local density increases and accre¬ 
tions are less likely to increase the spin value as much as 
major mergings at the same local density; meanwhile, a 
less massive halo tends to have a larger spin value at the 


same mass increase (Alogj^g-^)- regions, halos 

have the lowest value of Ac. From this description, we 
conclude that /td in lower density regions tends to be 
lower than that in higher density regions (or less-aligned 
infall directions), and accretion may favor a less biased 
infall direction lowering the characteristic spin. Also, 
smaller cluster halos tend to have more aligned mass in¬ 
falls (or a higher value of Ac) than larger cluster halos. 
The environmental dependence of Ac is more clearly ob¬ 
served in the case of higher Alog^Q M. 

The redshift dependence of Ac can be observed by com¬ 
paring Figure[5](0 < z < 0.2) with Figure[n](2 < z < 3). 
The dependence of Ac on the halo mass at a given local 
density is relatively negligible, especially for major merg¬ 
ing at high redshifts. Overall, Ac values at 2 < z < 3 are 
lower than those for 0 < z < 0.2 cases, especially for 
A logj^g M > 0. The environmental dependence becomes 
stronger at lower redshifts. 

We consider the frequencies of two infall scenarios 
(accretion or major merging) in Figure BOI Accretion 
(0 < AlogigM < 0.1) at low redshift (0 < z < 0.2, 
left panels) is more frequent than at the higher redshift 
(2 < z < 3, right panels). Moreover, halos in less dense 
regions at low redshifts are likely to have less frequent 
major-merging events; this environmental dependence is 
not observed at high redshifts. Therefore, we conclude 
that major-merging events subside with time, and this 
tendency becomes stronger in lower density regions. This 
may be due to the accelerating expansion of the universe 
at lower redshifts. The acceleration may reduce the num¬ 
ber of major-merging events; therefore, accretion would 
be a dominant source of mass increase in recent epochs. 
Lower-density regions, especially voids, may suffer more 
seriously from such acceleration. 

4.4. Random Walk and Correlated Mass Infall 

We need to clarify two confusing terms: random walk 
and correlated (or biased) infall. In this study, a ran¬ 
dom walk means that the angular momentum change at 
a given merging step is purely stochastic, obeying the dis¬ 
tribution of D. However, the D distribution (or P{D)) 
has, in most cases, has a hnite mean value: a positive 
bias {fiD >0) or negative bias (/td < 0). The correlated 
infall direction indicates a preferential direction for mass 
infall with respect to the halo rotation axis (or its spin 
axis; see Appendix |B] for details) and does not conflict 
with the random-walk hypothesis. 

Since /td is finite (and mostly positive), the orienta¬ 
tions of the orbital angular momentum of infalling matter 
are not perfectly random but are somewhat correlated or 
anti-correlated with the halo rotational axis (or the di¬ 
rection of the angular momentum vector). However, it 
should be noted that as the halo spin approaches zero, 
even a small amount of mass infall (or accretion) may sig¬ 
nificantly increase the halo spin value {\L\ :s> \J\ ~ 0), 
which may explain the sharp increase in /td as the sample 
spin approaches zero (Fig. [T]). 

An illustration of a spin random walk is a staggering 
walk taken by a drunken man. Consider a half-pipe with 
a U-shaped cross section and, then, a group of drunken 
men take a walk lurching forward along the bottom of 
the half-pipe without mutual collisions. Here, the shape 
of the cross section is not fixed, and the pipe is not like 
P{D), which varies with sample parameters. Therefore, 
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Fig. 8 .— Dependence of Ac on merging mass (Alogj^o-^) 
tween 0 < z < 0.2 for three mass samples: 0.3 < M 12 < 0.5 
(bottom panel), 1 < M 12 < 2 (middle), and 7 < M 12 < 10 (top). 
Symbols with different colors are used to distinguish the effect of 
local environment. 
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Fig. 10.— Distribution of dlogj^Q M for halo samples of 0.6 < 
Mi 2 < 0.8 (bottom panels) and 6 < M 12 < 10 (top) at 0 < < 0.2 

(left panels) and 2 < 2 < 3 (right). 


the distribution of walker positions gradually changes to 
adjust itself to the cross-sectional shape of the half-pipe, 
and the majority of walkers are distributed near the bot¬ 
tom of the pipe (Ac in the halo spin walk). In the next 
section, we investigate the results of spin random walks 
and the effects of various parameters on the shape of spin 
distributions. 

5. RANDOM EVOLUTION OF HALO SPIN 

In this section, we will show how well the stochastic 
model reproduces distributions of halo spin measured 
from V-body simulations. The change in angular mo¬ 
mentum is expected to be a random process with a dis¬ 
tribution that has a bimodal Gaussian form. We assume 
that the spin change in each time step is ind ependent of 
the p revious history (a Markov process; see iBond et ^ 
II 99 II for an example applied to halo mass function) and 
is completely determined by the current physical status. 
However, it is possible that subsequent infall events can 
be correlated with previous events. For example, con¬ 
sider successive mass infalls through a local filamentary 
structure. In this case, the random-walk hypothesis is 
not satisfied. We will discuss possible correlated infalls 
in the latter part of this section. 

5.1. Stochastic Model of Spin Evolution 

First, we choose a mass history from simulated main 
merger trees. The initial spin value is set to Aq = 0.01. 
At each merging step, we identify P{D) of a subsample 
with parameter ranges enclosing the current halo mass, 
spin, mass change, local density, and redshift. Then, 
from the identified P{D) we randomly generate a value 
of T^rand- From equation (l9|) with Hrand, we calculate the 
change in spin (AA) and update the spin value. In the 
next merging step, we iterate the same process but with 
the updated spin value. The resulting spin value (A/) at 
the final redshift is obtained as 

N 

= ^0 P AA^. 

i=l 


Fig. 9.— Same as Fig. [8] but for redshift samples of 2 < 2 < 3. 


(16) 

























9 


5.2. Effect of Halo Mass on Spin Distribution 

Now, we obtain the spin distribution produced by the 
stochastic random walk and compare it with that of 
7V-body simulations. Figure [TT] shows the results with 
randomly generated spin distributions (blue histogram) 
against 7V-body simulated distributions (red histogram) 
in four mass samples. The randomly generated spin dis¬ 
tributions at 2 = 0 are similar to the simulated distribu¬ 
tions with a slight shift to higher A, and this difference is 
observed in every mass sample, although the difference 
is negligible at higher redshift. 

In Table |4] in Appendix m the log-normal fitting re¬ 
sults for randomly generated and simulated spin distri¬ 
butions are given with a mean value, standard deviation, 
and x^/degrees of freedom (dof). This difference is more 
clearly shown in Figure [121 where the blue and red his¬ 
tograms are the distributions of the random and simu¬ 
lated spins at z = 0, respectively. The green curve is a 
Gaussian fit to the simulated spin distribution. Before 
jumping to a detailed comparison, it would be valuable 
to note that iV-body-simulated halos overpopulate at the 
high-A tail (the same feature is also observed for the gen¬ 
erated spin distributions) compared to the log-normal 
fitting ( solid curve), which has also been f ound in other 
studies (|Bett et al.ll2007l : fShaw et al.ll2006ll , unless a cor¬ 
rection is applied to the halo samples (iDavis et al.ll2011l : 
IShaw et al.(l2006f) . In Figure [T31 we show the difference 
= {Pr — Ps)/Ps) between the simulated (Pg) and 
randomly generated (Pr) spin distributions. The differ¬ 
ences increase as redshift decreases, while deviations in 
the randomly generated spin distributions at z = 0.5 and 
2 are confined within about 10-15% error. 


that a randomly generated spin (blue histogram) tends to 
substantially overestimate the halo spin of the A^-body- 
simulation in lower density regions. 

For the group and cluster samples, the difference 
in spin distributions between the simulation and the 
stochastic model is relatively small. However, with time 
the difference increases with time and the group sample 
begins to show such a discrepancy at z = 0. 

Since a massive halo is likely to be located in a knot (in¬ 
terconnection) of local filamentary structures, the proba¬ 
bility of successive major mergers along a single filament 
is low (the next major merger may occur with another 
filamentary structure). Therefore, it would be interest¬ 
ing to determine whether the randomly generated halo 
spins of more massive samples have a distribution simi¬ 
lar to that of the A^-body simulation. In Figure [T71 we 
show the generated spin distributions of a massive sam¬ 
ple of 6 < Mi 2 < 10 in different local environments. 
No substantial difference is shown even in void regions. 
Therefore, it is quite reasonable to believe that massive 
halos tend to show random evolution in spin. 

In these figures, the halo spin distribution depends on 
the local environment; denser regions cause halos to ro¬ 
tate faster, which is also confirmed in Figure [5] (Ac is 
higher in a denser region). Also, most of the spin sam¬ 
ples show log-normal distributions, except for the cluster 
samples at z = 0, where the distribution has a slightly 
fat tail at larger A. This may be due to unrelaxed halos 
in cluster regions where violent mergers are so frequent 
that cluster halos do not have sufficient time to transfer 
angular momentum to outside and finally to be virialized. 

6 . MARKOVIAN AND NON-MARKOVIAN PROCESSES 


5.3. Environmental Effect on Spin Distribution 

If a spin walk is well described by the Markov (or 
random) process and a sufficient number of trials are 
performed to suppress the Poisson error, then the N- 
body-simulated and randomly generated spin distribu¬ 
tions should be basically equal to each other. Therefore, 
a subtle failure in recovering the simulated spin distribu¬ 
tion at lower redshifts, as seen in Figure [11] may imply 
a possible limitation of the stochastic model. 

A possible explanation for the mismatch of spin distri¬ 
butions between the simulation and the stochastic model 
is correlated mass merging, which means a subsequent 
mass infall changes the angular momentum of a halo in 
the same direction. Then, a set of lArandS in the suc¬ 
cessive infalls may not be random but correlated. This 
correlated mass merging could be a function of local en¬ 
vironment. In void regions, halos are likely to have fewer 
major mergers, while those in cluster regions may experi¬ 
ence more frequent heavy mass infalls (see Fig. [TOj). Be¬ 
cause it is believed that major mergers frequently occur 
when clumps of matter fall along filamentary structures, 
the number of filaments connected to a halo determines 
whether or not the spin walk is random. In void regions, 
filamentary structures are relatively rare and void halos 
tend to have a small number of filaments around them, 
which means that subsequent major mergers may occur 
along a smaller number of filamentary structures. 

Figures [HHTBl show the environmental effect on spin 
distribution at three redshifts. In Table [5| of Appendix [F] 
the log-normal fitting results are listed. We clearly see 


Throughout this paper, the stochastic model we 
adopted assumes that the spin walk can be described 
by the Markov process, which requires that a random 
walk be memoryless or independent of previous states. 
Therefore, it is important to determine whether the spin 
walk observed in A^-body simulations is really random 
(Markovian) or not (non-Markovian). 


6.1. Autocorrelation between Neighboring Mass Events 

One of the easiest ways to verify the correlation of two 
separate events in a halo merging history is to measure 
the autocorrelation (or serial correlation) between a pair 
of events in a main merger tree line. The first-order 
autocorrelation of two mass events of i and j is 


l{i,j) = WiWj 


Eli - pii\ (Dj - fij 


(17) 


where re is a weight, D is randomly generated by the 
stochastic model or provided by simulations, pi is the 
mean, and cr is the standard deviation of the prob¬ 
ability distribution of D. In this analysis, we adopt 
Wi = |Alog]^gMi| and define the lag-fc autocorrelation 


as 


Gfc(z) 


1 

Wpair 


^{i,j\>z), 


(18) 


where Wpair = WiWj and the summation is performed 
over all possible lag-/c pairs of events in a halo history to 
redshift z. Only a relative comparison is possible because 
the time interval between a lag-fc pair changes. 
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Fig. 11.— Spin distributions of various halo mass samples at 2 : = 0, 0.5, and 2 (from left panels). The blue and red histograms are the 
randomly generated and A^-body-simulated spin distributions, respectively. 


Using 100 random realizations of P(Cfc), we statisti¬ 
cally measured the deviation of the autocorrelations of 
iV-body data from random expectations. The proba¬ 
bility distribution of a lag-one autocorrelation is shown 
in Figure [15] for halo samples of 1 < M 12 < 2 in four 
local-density regions (colored symbols). The brown his¬ 
togram is a mean expectation measured from the 100 
random realizations. Even though the probability distri¬ 
bution of D depends on the local density, the autocor¬ 
relation seems to be indistinguishable between different 


local-density samples and even from the random model 
in this plot. 

We also showed deviations between simulated corre¬ 
lations and the random model with ^P(Cfc) = {Ps{Ck) - 
P<r>{Ck))/P<r>{Ck), where Ps{Ck) is the simulated au¬ 
tocorrelation and P<r> (Ck) is the average of the 100 ran¬ 
domly generated autocorrelations. Figure [H] shows the 
results for a sample of 1 < M 12 < 2. The shaded re¬ 
gion depicts the Icr distributions of the random model, 
and the symbols are measured from 7V-body simulations. 
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Fig. 12.— Spin distributions of the mass sample of 1 < M 12 < 2 
at 2 : = 0. Blue and red coded histograms are the randomly gen¬ 
erated and A^-body spin distributions, respectively, and the green 
curve shows the log-normal fit to the A^-body histogram. 



Fig. 13.— Deviations in spin distribution for three mass samples 
measured in the spin interval of 0.01 < A < 0.1. The circles, boxes, 
and triangles are measured at 2 : = 0, 0.5, and 2, respectively. We 
removed the irrelevant regions of A < 0.01 and A > 0.1 as the 
number of simulated halos in those regions is very small. 

It is easy to see the negative autocorrelations of Ci at 
z = 2 and the positive correlations of 02 at 2 ; = 0 with 
greater than 2a confidence limits. We further investi¬ 
gated autocorrelations from C 4 to Cq but found nothing 
substantially deviating from the random model. On the 
other hand, we could not detect any possible autocorre¬ 
lations in the massive sample (6 < M 12 < 10)). Most of 
the simulated autocorrelations are contained within Icr 
random expectations (see Fig. HOD. 

Detecting substantial autocorrelations of Ci at z = 2 
and (72 at z = 0 may provide a tangible clue for void 



Fig. 14.— Spin distributions in various local environments at 
2 : = 2. Counterclockwise from the bottom left panel are the spin 
distributions of A^-body (red histogram) and random-generated 
(blue) samples of local densities of Apio < 0.7, 0.7 < Apio < 2, 
2 < Apio < 10, and 10 < Apio < 100. The green solid curve in 
each panel is a log-normal fit to the A^-body spin distribution. 



Fig. 15.— Same as Fig. 1 141 but at 2 : = 1. 

discrepancies between the simulated and generated spin 
distributions. Correlated infall events may cause void 
halos to have smaller spin values, but we are unable to 
detect any environmental dependence in the correlation 
hgure. However, considering the frequent major merging 
of cluster halos, an autocorrelation could easily be coun¬ 
terbalanced by other major-merging events. However, 
the void halos have little possibility of another major 
merging. 

6 .2. Origin of the Log-normal Distribution of Halo Spin 

In this subsection, we demonstrate that the log-normal 
distribution is a natural consequence of the stochastic dif¬ 
ferential equation in the stochastic of halo spin. In many 
studies, simulated spin distributions in various mass sam¬ 
ples consistently show substantial deviations from the 
log-normal form in the high-A tail. However, no seri¬ 
ous deviations are found as long as a mass sample is 
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Fig. 16.— Same as Fig. 1141 but at z = 0. 



Fig. 17.— Spin distributions of a more massive sample of 6 < 
Mi 2 < 10 at z = 0. The green solid curve is a log-normal fit to the 
corresponding Al-body distribution. 


further divided into local-density subsainples. As shown 
in Figures [1^161 the log-normal fitting functions (green 
curves) seem to describe well the simulated spin distri¬ 
butions (red histograms) in various local environments, 
except for cluster regions at z = 0. 

Now, consider a one-factor model of the Markov pro- 
cess or the model of Geometric Brownian Motion (GBM; 
IRossI[2 M01 . which has a stochastic differential equation 
of 


dlogio A(r) 
dr 


dWr 

dr 


(19) 


where 6 describes the long-term drift of the system, fjc 
is a constant, and Wr is a kind of normally distributed 
Wiener process or Wr ^ A/*(0,r) where A/*(0,r) denotes 
the normal distribution with a zero mean value and a 
variance of r. 

From the properties of the Wiener process such as 
0 ^t 2 — ^Ti) A/*( 0 ,T 2 — Ti) for r 2 > Ti, we obtain 



c, 


Fig. 18.— Example of the probability distribution of Cl fEg. llSfl 
in a sample of 1 < M 12 < 2 at 2 : = 0. The symbols are distribu¬ 
tions measured with the spin history of the A^-body halo sam¬ 
ple with Apio < 0.7 (green), 0.7 < Apio < 2 (magenta), and 
2 < Apio < 10 (blue). The histogram shows the mean distribu¬ 
tion of 100 random generations in the sample of Apio < 0.7. 


the probability distribution of dW/dr as 


dWr 

dr 


Wr+dr - Wr 
dr 


^(0, dr) 
dr 


= V(0,l/dT), (20) 


where the statistical relation of = 

J\f/a^) is used. Thus, the probability distribution 
of dWr /dr is a normal distribution with a zero mean and 
a variance of l/dr. The Wiener process has the following 
properties: 


dWA-A/'(0,l)vT^ (21) 

dW'^ ~ vTV(0, l)dT oc dr. (22) 

According to Ito’s formula (|Movellan. J. R.ll2011[l . we 
may linearize the change of log A as 


dlog A 


dlogA^^ , 


= acdWr + 



(23) 

(24) 


where we used the relation dA = \{9dT -F acdWr) (from 
Eq. [T9|) and dWr^ = dr (Eq.[2pl), and ignored the higher- 
order terms of dr. By integrating the above equation, we 
finally obtain 

Now, one can understand that the logarithm of the 
stochastic process A follows a Wiener process with a stan¬ 
dard devi ation of ar and a corrected long-term drift of 
Q — all2 (|0ksendall (j200nH for details). If W-r does not 
have a normal distribution, then we may not obtain a 
logarithm of the randomly distributing quantity. Eor ex¬ 
ample, if Wr has a probability proportional to —r or —r^, 
then the resulting distribution would be exponential or 
Gaussian, respectively. 
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Fig. 19.— Randomness check by the measurement of autocorrelations, Ci (left panels), C 2 (middle), and C 3 (right) in the sample of 
1 < Mi 2 < 2 at 2 = 0 (bottom panels) and 2 (top). The three gray levels in the regions are the 1-a distributions of 100 random generations 
for Apio < 0.7 (outer gray region), 0.7 < Apio < 2 (middle), and 2 < Apio < 10 (innermost) with corresponding A^-body estimations 
plotted with green, magenta, and blue circles, respectively. 



Fig. 20.— Same as Fig. I19l but for the massive samples of 6 < M 12 < 10. Due to the smaller sample size, the 1-cr scatter of the random 
sample is larger than that of the sample of 1 < M 12 < 2 . 
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We now attempt to combine the right two right-hand 
terms in equation (HI into V as 


dlogio A(t) 
dr 


(26) 


where V ^ dr]. After changing variables as 

r —>■ log^Q M and 1) ^ D — Dc, we can find an analogy 
between the GBM and the angular momentum change, 
D (see Eq. O. Although the overall distribution of D 
is better fit by the bimodal Gaussian, the distribution 
around a peak can be well modeled by a single Gaussian 
peak because the two Gaussian components have nearly 
the same center in most samples, and we may simply 
approximate P{D) around the peak as a normal distri¬ 
bution. The minor Gaussian component with a broad 
width is mostly needed to describe the fat tails on both 
sides. Gonsequently, the distribution of log^o A follows 
the normal process, which means that the spin distribu¬ 
tion of A should accordingly be log-normal. If the centers 
of two Gaussian components are well separated, or the 
minor Gaussian component is comparable to the major 
one, then the resulting distribution of log^g A may show 
a considerable deviation from the normal distribution. 

Here, we need to address the validity of the application 
of the GBM model to spin evolution before jumping to 
a conclusion, and thus we test whether the halo spin 
change satisfies two prerequisites of the GBM model: 


• 9 (or (Id) should depend only on r (or log^o-M); 
and 



• the standard deviation of V (or ctd) should be pro¬ 
portional to (or (A logio 

To qualitatively assess the first condition. Figures [31 HI 
|5l and [6] show that there is no significant change in 
with different redshifts, local environments, or infall mass 
ratios, except for major-merging events, in which case a 
substantially larger /iu is obtained. Consequently, one 
can easily see the mass dependence in Figure |7l 
The second condition might be related to the tempo¬ 
ral resolution affecting the mass infall ratio. Now, we 
investigate the following relation: 

crjji = (AlogioM)^/2^-\ (27) 

where CTc is constant with respect to A log^o M as long 
as the second requirement is satisfied. In order to obtain 
a log-normal distribution of A, should be linearly 
scaled to the change in halo mass. As shown previously, 
the spin distribution has a mass and environmental de¬ 
pendences, indicating that the distribution is not station¬ 
ary but evolves with the halo mass. However, halo envi¬ 
ronments do not change substantially over cosmic time 
since a halo does not generally move a great distance (< 
a few Mpc). Hence, it is sufficient to examine the relation 
between cr~jy and Alogj^g M in the same environment. 

Figure |2T] shows the change in over a range of 
infall mass ratios at z = 0. The measured relation sub¬ 
stantially deviates from a single power-law scaling but 
we can divide the range of A logj^g M into three regions: 
Alog^gM < 0.03 (quiet accretion), 0.03 < Alog^o-M < 
O.I (strong accretion), and A log^g M > O.I (major merg¬ 
ing). Each region, except for the strong accretion mode, 


Fig. 21. — Changes in and infall mass ratio, Alog^gM 
(bottom), at 0 < x: < 0.2 for different halo masses, 0.3 < M\2 < 
0.5 (bottom), 1 < Mi 2 < 2 (middle), and 5 < M 12 < 7 {top 
panel). The dotted lines show the linear relation between and 
d logic 

has its own power-law relation of oc (Alogig M)^/^. 
The major-merging mode and quiet accretion mode have 
different variances in that the major-merging mode has 
a higher value of cTc- This implies that the major merger 
mode has a relatively larger stochastic walk size than the 
quiet accretion (see Eq. dH). The transition (or strong 
accretion) mode (0.03 < A logig M < O.I) does not seem 
to satisfy the assumptions of the GBM model. Now we 
conclude that even though there exist limitations in ap¬ 
plying the GBM model to spin evolution, the spin walk is 
well described by the GBM model and, consequently, the 
resulting spin distribution follows the log-normal func¬ 
tion. 

7. CONCLUSIONS 

We concluded that the spin distribution is well de¬ 
scribed by a random walk of angular momentum. Halos 
in lower density regions in recent epochs present some 
failures. This may be due to correlated major-merging 
events (see further discussions in Appendix iBl and (Cl). 
The log-normality of the spin distribution is found to be 
a consequence of the stochastic random walk, and most 
density samples show log-normal distributions of the halo 
spin. Only a small departure is found in the cluster sam¬ 
ples at z = 0, mainly because a fraction of halos are not 
yet relaxed. We found that the angular momentum of 
a halo is likely to increase after merging, and the spin 
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distribution depends on the sample mass. If the sam¬ 
ple mass is higher, then the average spin value decreases. 
The simulated log-normal distribution of the spin is well 
recovered by the stochastic model with bimodal Gaussian 
distributions of angular momentum change. 

In the standard ACDM model, the universe has un¬ 
dergone a recent accelerating expansion, making it dif¬ 
ficult to pull nearby material into a local shallow gravi¬ 
tational center. Compared with the flat Vim = 1 model, 
the ACDM model has prominent filamentary structures 
lacking many small and faint structures that could be 
easily destroyed by the recent cosmic acceleration. This 
effect is even stronger in void regions where the mass 
evolution was already halted at higher redshifts. There¬ 
fore, the cosmological effects may leave clear evidence in 
the spin evolution of void halos and mean-field halos at 
lower redshifts (0.7 < Apio <2; see the top-left panel of 
Fig. ITCl) . 

Because of the bimodal Gaussian shape of the angular 
momentum change in the stochastic model, one may raise 
the question of the possible existence of another hidden 
parameter for subsampling merging events to measure 
P{D). By introducing another parameter, a simpler dis¬ 
tribution shape may be obtained. A bimodal Gaussian 
shape could be obtained through a possible correlation 
among spin walks, which seems to be signihcant in void 
regions and in mean helds. This may explain the de¬ 
viations in halo spin distributions from the log-normal 
distribution. 

In this study, we directly adopt the simulated mass 
evolution for the halo mass change (Alogj^gM). There 
are several mod els for the mass merging based on the 
EPS fo rmalism (iPress &: Schechted (1197^ for the classi- 


nce^.^.Col€ 


cal an d [Bower] (|1 991li: Bo nd et al.l dlQQlfl : [ 

(1993); Kaiiffmann fc White! (1199311: iMo fc White 
199(B: ISheth fc LenrsonI (1199911: Sheth et al.1 


Sheth fc Tormeni (1200^: I van den Boschl 

Hiotelis fc Del Popolol ( 20061) : Zhang et al.l 

iNeistein fc Dekeil 


Moreno et al. 


( 20011 ) 

( 2002 

(200g 

(200g 


Parkinson et al.l (j2008[i : l.liang fc van den Boschl l)201 


for variant EPS models), which can also provide us 
the probability of merging phases (major merger or 
accretion) between time steps. Since our stochastic 
model may discriminate between the effects of accretion 
and major merging in spin distributions, we may 
check whether the variant EPS models may produce 
proper probability distributions of AM/M compared to 
simulated ones for a given time step. In Appendix El 
we show the roles of major merging and accretion in 
shaping the spin distribution in detail by adopting a toy 
model with fixed the mass-merging ratio and how much 
the mass-merging history generated by the classical 
EPS model predicts different spin distributions from the 
A-body simulated ones. 
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Fig. 22.— Effect of 0 ( 2 :) on the spin distribution at z = 1. 
The histograms colored in blue and magenta are obtained with 
and without a{z) in Eq. [S] respectively. As references, the red 
histogram and a green curve are the A-body spin distribution and 
its fitting function, respectively. 
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APPENDIX 

EFFECT OF a{Z) ON SPIN DISTRIBUTION 

a(z) deserves more investigation because its effect on 
spin distribution is not negligible. The cosmological con¬ 
tributions to log-spin change with a mass increase in 
Alog^o M are —0.008 < a( 2 ;)Alog 2 o M < 0 in most red- 
shift intervals of merging data. So, even though its con¬ 
tribution is neglected, spin error is less than 1%. But, 
the error accumulates over the number of merging steps 
and may substantially change the spin distribution. 

Figure [H] shows how the spin distribution changes 
when the a(z) term is neglected in the random walk 
process. The generated spin distribution is substantially 
shifted to higher A, which demonstrates the sensitivity of 
spin walk to such a very small value or possible error in 
fitting P{D) and how large the sample size should be to 
suppress this kind of a noise effect. However, this does 
not guarantee that we may determine the cosmological 
model from the measured spin distribution because the 
measured P{D) in this analysis may differ in different 
cosmological models. 

ANGLE ALIGNMENT BETWEEN ROTATIONAL AND 
ORBITAL ANGULAR MOMENTA 

It would be interesting to examine the angle alignment 
between the rotation (J) of a halo and the orbital an¬ 
gular momentum (X) of infalling matter. If they are 
positively or negatively aligned, then the rotational an¬ 
gular momentum of a halo is expected to increase or de¬ 
crease, correspondingly. Therefore, this analysis is com¬ 
plementary to the work on the change in angular mo¬ 
mentum [D). In this section, we study two infall modes: 
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accretion (0 < Alogj^p-^ < 0-1) major merging 
(Alog^o^f > 0.1). The angle between two vectors is 
measured by 

^ ( 1 ^) ■ 

Figures [23l [24l and [25] show the probability distribution 
of 9a in halo samples of 0.03 < A < 0.05 (low-spi), 
0.01 < A < 0.03 (average-spin), and 0.05 < A < 0.3 
( high-spin), respectively, for two mass-merging modes 
(the accretion in the left and major merging in the right 
panels). In the accretion mode, the low-spin halos (Fig. 
[Ml) are likely to have slightly negative correlations with 
9a (anti-aligned), while average-spin halos tend to have a 
slightly positive alignment, and high-spin halos (Figs. (Ml 
and [25|l have a strong positive alignment. On the other 
hand, in major-merging mode, the low-spin, average- 
spin, and high-spin samples show negative, random, and 
positive correlations, respectively. 

If two angular momentum vectors are anti-aligned, 
then a halo’s rotational angular momentum (J) is re¬ 
duced and, consequently, its spin value is lowered. There¬ 
fore, a halo with an average spin value is likely to have 
a slightly positive correlation to offset the negative drag 
term (—5/3), while a high-spin or low-spin halo tends 
to have positively or negatively aligned mass infalls, re¬ 
spectively. Then, the alignment of two vectors and the 
value of a halo spin are closely correlated with each other, 
producing a high spin value, and vice versa. 


ANGLE CORRELATIONS BETWEEN MASS EVENTS 


This section discusses testing to determine whether or 
not the directions of the orbital angular momenta of in¬ 
falling matter in two events are correlated. The angle 
between the orbital angular momenta in the ith and jth 
events is defined as 


9ij = cos 




(Cl) 


and we apply a weighting of = A logj^g M placing 
more weight on major-merging pairs. We only consider 
j = i + 1 cases. Figures (Ml- (Ml show the directional cor¬ 
relations of the slowly, mildly, and fastly rotating halos, 
respectively. In the figures, there is no explicit depen¬ 
dency on environment. Slow-rotating halos show a neg¬ 
ative correlation (anti-aligned) with regard to the direc¬ 
tion of L between two successive events, which explains 
why they are slowly rotating. However, some interesting 
features can be found in the mildly and fastly rotating 
samples, which develop a strong anti-correlation at lower 
redshift (z < 1; bottom-left panel of Figs. iTTl and (Ml). 
This may curb the development of spin as a whole and 
may explain the reason why void halos seem to have a 
non-Markovian walk of D. Halos in denser regions are 
expected to experience frequent major mergers. Thus, 
even though there is anti-alignment between two steps, 
the many other major merger events violently erase the 
memory of previous anti-correlations. Therefore, the su¬ 
perpositions of multiple anti-correlations are completely 
mixed and may result in a random-like process. How¬ 
ever, void halos tend to experience relatively few major 
mergers, and the effect of anti-alignment in a rare pair of 
merging events may survive for a longer time. This may 



0 50 100 150 0 50 100 150 


0A 0. 


(a) accretion 



0 50 100 150 0 50 100 150 


0A 0. 


(b) merging 


Fig. 23.— Angle alignment between J and L of accretion (a) in 
major-merging events and (b) in halo samples of 0.03 < A < 0.05 
and 1 < Mi 2 < 2. Counter-clockwise from the lower left panel, 
sample redshifts are 2 : = 0, 1, 2, 3, 4, and 5. Each color-coded 
curve is measured from the local-density sample of void (blue), 
mean-field (red), or group (green) regions. The magenta sine curve 
(dotted line) marks the random orientation. 
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Fig. 24.— Same as Fig. 1231 but for slowly rotating samples of 
0.01 < A < 0.03. 



Fig. 25.— Same as Fig. 1231 but for fast-rotating samples of 
0.05 < A < 0.3. 


explain why void halos are not properly described by the 
stochastic random walk of D. We discuss the number of 
major mergers a halo may experience during its life time 
in Appendix lEl 


AUTOCORRELATIONS BETWEEN MA,JOR MERGERS 

We estimate the autocorrelation between a pair of 
nearest major merger events to determine whether sub¬ 
stantial amounts of correlation could be detected. In this 
case, we only consider autocorrelations between major 
merger events. The deviations of autocorrelations be- 
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Fig. 26.— Angle correlations (^ 12 ) of orbital angular momenta 
for subsequent mergings in slow-rotating samples with 0.01 < A < 
0.03. Each color-coded curve is measured from the local-density 
sample of void (blue), mean-field (red), or group (green) halos. The 
magenta sine curve (dashed line) marks the average expectation of 
random orientations. 



Fig. 27.— Same as Fig. 1261 but for the mildly rotating sample 
with 0.03 < A < 0.05. 



Fig. 28.— Same as Fig. [26] but for the fastly rotating sample 
with 0.05 < A < 0.3. 


ers, which may explain the deviations of the stochastic 
model prediction from A^-body spin distributions. In Fig¬ 
ure [301 we show the number of major mergers (Vmm) per 
halo with masses 1 < M 12 < 2 (a) or 6 < M 12 < 10 (b). 
In the sample of 1 < M 12 < 2, the number difference 
between halos in different environments is not severe at 
higher redshift {z > 2), while, at lower redshifts, less 
dense samples clearly show fewer major mergers. How¬ 
ever, more massive samples (b) show nearly negligible 
deviations from each other, demonstrating noticeable dif¬ 
ferences only at z = 0. The peak number of the lYmm 
distribution is four, while it is three at higher redshifts 
(z > 4). Therefore, we conclude that the correlated ran¬ 
dom walk seen in void regions may be a joint effect of the 
major-to-major merger autocorrelation and the smaller 
number of major mergers in less dense regions. It is in¬ 
teresting to find that a significant fraction of halos in 
cluster regions do not experience a major merger event 
during their lifetime (5% of halos with 1 < M 12 < 2 and 
1 % of halos with 6 < M 12 < 10). Also, the fraction 
of cluster halos needed to prohibit a major merger event 
after z = 1 is significantly higher than that of halos in 
lower density regions. 


tween major mergers (Cmm) are shown in Figure |29| for 
several redshifts. In this figure, major merger autocor¬ 
relations at high redshifts (z > 2) are equivalent to the 
random distributions, while weak positive correlations 
are detected (~ Icr) at lower redshifts. More massive 
samples do not show any positive or negative correla¬ 
tions, as seen in panel (b) of Figure |29l Therefore, in 
this plot, we are sure that there is no substantial auto¬ 
correlation between major mergers. 


STATISTICAL PROPERTIES IN THE SPIN 
DISTRIBUTION 

In this section, we present the means and standard 
deviations of spin distributions from V-body simulations 
and random generations together with fitting results to 
a log-normal function with least . See Tables [4] and 
m The mean value of the halo spin increases with local 
density and redshift. The difference between V-body and 
random simulations increases as local density decreases. 


NUMBER DISTRIBUTIONS OF MA.IOR MERGERS 

We determine the distribution of major merger fre¬ 
quency in each halo sample at z = 0 to confirm our 
argument that void halos experience fewer major merg- 


EFFECT OF MERGING MODEL ON SPIN DISTRIBUTION 

In this section, we demonstrate the effects of major 
merging and accretion on the spin distributions in a sim¬ 
ple toy merging model. The mass increase is parameter- 
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TABLE 4 

Spin Distributions for the Halo Mass Samples Used in Figures 




A-body Simulation 

Random Generation 

Halo Mass 

Kedshitt 

(logio A) 

“"loglO ^ 

X^/dof 

(logio A) 

“"logio I' 

X^/dof 


0.0 

-1.386 

0.336 

0.00118 

-1.342 

0.329 

0.00140 

0.5 < Mi2 < 0.8 

0.5 

-1.389 

0.319 

0.00076 

-1.374 

0.308 

0.00071 


2.0 

-1.434 

0.290 

0.00034 

-1.429 

0.289 

0.00037 


0.0 

-1.407 

0.323 

0.00068 

-1.358 

0.326 

0.00107 

1 < Mi2 < 2 

0.5 

-1.405 

0.311 

0.00037 

-1.381 

0.309 

0.00059 


2.0 

-1.450 

0.291 

0.00029 

-1.427 

0.294 

0.00031 


0.0 

-1.414 

0.313 

0.00038 

-1.366 

0.324 

0.00080 

4 < Mi2 < 6 

0.5 

-1.412 

0.304 

0.00028 

-1.386 

0.306 

0.00049 


2.0 

-1.474 

0.291 

0.00064 

-1.439 

0.296 

0.00062 


0.0 

-1.415 

0.309 

0.00027 

-1.367 

0.322 

0.00090 

7 < Mi2 < 10 

0.5 

-1.419 

0.302 

0.00028 

-1.390 

0.306 

0.00061 


2.0 

-1.488 

0.292 

0.00086 

-1.446 

0.304 

0.00110 


Note. — The mean and standard deviation of the halo spin from N’-body simulation and 
random generation. Also, we list the least obtained during the log-normal fitting. 


TABLE 5 , 

Spin Distributions for the Halo Samples Used in Figures HaITtI 



A-body Simulation 

Random Generation 

Kedshitt Environment 

(logic A) “login A xVdof 

(logic A) “login A X^/dof 


1 < Mi2 < 2 



void 

-1.539 

0.277 

0.00108 

-1.440 

0.297 

0.00083 

z = 0 

mean-field 

-1.486 

0.284 

0.00041 

-1.424 

0.297 

0.00057 


group 

-1.410 

0.311 

0.00050 

-1.357 

0.327 

0.00083 


cluster 

-1.287 

0.375 

0.00238 

-1.287 

0.374 

0.00130 


void 

-1.544 

0.273 

0.00125 

-1.432 

0.294 

0.00059 

Z = 1 

mean-field 

-1.489 

0.279 

0.00053 

-1.435 

0.291 

0.00055 


group 

-1.412 

0.293 

0.00028 

-1.391 

0.299 

0.00050 


cluster 

-1.317 

0.320 

0.00058 

-1.332 

0.320 

0.00065 


void 

-1.564 

0.272 

0.00168 

-1.457 

0.290 

0.00081 

z = 2 

mean-field 

-1.522 

0.277 

0.00069 

-1.461 

0.289 

0.00043 


group 

-1.454 

0.285 

0.00036 

-1.424 

0.294 

0.00038 


cluster 

-1.369 

0.298 

0.00026 

-1.367 

0.303 

0.00043 


6 < Mi2 < 10 



void 

-1.574 

0.268 

0.01237 

-1.529 

0.319 

0.01626 

z = 0 

mean-field 

-1.516 

0.285 

0.00112 

-1.459 

0.318 

0.00145 


group 

-1.423 

0.300 

0.00036 

-1.378 

0.315 

0.00073 


cluster 

-1.335 

0.328 

0.00072 

-1.290 

0.347 

0.00101 


void 

-1.635 

0.274 

0.03620 

-1.585 

0.305 

0.02587 

2 = 1 

mean-field 

-1.549 

0.280 

0.00339 

-1.504 

0.296 

0.00257 


group 

-1.447 

0.288 

0.00074 

-1.421 

0.294 

0.00067 


cluster 

-1.360 

0.298 

0.00062 

-1.335 

0.301 

0.00090 


void 

-1.679 

0.283 

0.10107 

-1.589 

0.349 

0.10899 

2 = 2 

mean-field 

-1.613 

0.270 

0.01070 

-1.537 

0.317 

0.00637 


group 

-1.520 

0.281 

0.00169 

-1.472 

0.302 

0.00155 


cluster 

-1.435 

0.293 

0.00154 

-1.384 

0.306 

0.00190 


Note. — Same as Table |4] 

ized as 


where /3 is set constant over the halo mass evolution. The 
number of time steps is, therefore, determined from the 
input value of /3 for given initial and final halo masses. In 
the top right panel of Figure[3Tl we observe that a higher 
value of P makes the spin distribution move to a higher 
A. However, the actual mass growth is a spectrum of 
AM/M with a certain probability distribution and this 
shapes the spin distributions found in simulations. 

Next, we test whether the mass accretion histories de¬ 


rived from the classical extended Press-Schechter (EPS) 
formalism may give a probability distribution of AM/M 
comparable to the simulated ones. We generate merg¬ 
ing trees according to the EPS merging mode l of the 
spherical top-hat collapse (jLacev fc Colei 1199,111 . The 
conditional probability of a progenitor to have a den- 
sity contrast of Srjt^) at the smoothing scale of Mi is 
(iJiang fc van den Boschl[2014ll 

/sc(5'i, WilS'o, Wo) = --J= ’ (^2) 

where AS = Si — Sq, Si = cr'^{Mi), Aw = wi — wq, and 
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Fig. 29.— Autocorrelations between the major mergers at red- 
shifts 2: = 0, 1, 2, 3, 4, and 5 in the samples of 1 < M12 < 2 (a) 
and 6 < M12 < 10 (b). 
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Fig. 30.— Number distributions of major mergers in the history 
of a halo until 2: = 0, 1, 2, 3, 4, and 5 measured from the samples 
of 1 < Afi2 < 2 (a) and 6 < M12 < 10 (b). 


LUi = 6c(ti) for fo > ti- Using this equation, we gen¬ 
erate the 100,000 random mass-merging trees and follow 
the stochastic spin changes for several halo mass samples 
(see Fig. [3T|) at ^ = 0. In the figure, the EPS model pre¬ 
dicts that the spin distributions will shift to higher values 
of A, but this tendency is smaller for larger mass samples. 
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Fig. 31. — Spin di stributions from the EPS model of the spheri¬ 
cal top-hat collapse JLacev fc Coleiri99,H l: top-left and bottom row 
panels) and from the toy model with various /3’s (top-right panel) 
at 2 = 0. 


This means that the spectrum of merger rates from the 
least accretion to most-massive major merger modeled by 
the classical EPS model may be different from 7V-body 
results. Thus, the results suggest that it would be inter¬ 
esting to test whether other variant, and certainly more 
advanced, EPS models that have been well known to 
provide halo mass functions consistent with simulations 
may also produce spin distributions similar to those from 
A^-body simulation by adopting this stochastic modeling. 
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